Finite-Size Scaling at the Jamming Transition 
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We present an analysis of finite-size effects in jammed packings of N soft, frictionless spheres at 
zero temperature. There is a i correction to the discrete jump in the contact number at the tran- 
sition so that jammed packings exist only above isostaticity. As a result, the canonical power-law 
scalings of the contact number and elastic moduli break down at low pressure. These quantities 
exhibit scaling collapse with a non-trivial scaling function, demonstrating that the jamming transi- 
tion can be considered a phase transition. Scaling is achieved as a function of N in both 2 and 3 
dimensions, indicating an upper critical dimension of 2. 



Numerical simulations of particulate systems are un- 
avoidably limited to a finite number of particles. It 
has long been recognized in the context of phase tran- 
sitions that this limitation can be exploited [1] to yield 
insight into the nature of the transition. In the context 
of the zero-temperature jamming transition of friction- 
less sphere packings [2], a finite-size analysis should be 
particularly valuable because the transition appears to 
be a rare example of a random first-order transition in 
finite dimensions, characterized by a discontinuous jump 
in the contact number (i.e. the average number of inter- 
acting neighbors) and power-law scaling [3, 4] as well as 
diverging length scales [5-7]. 

In this paper, we establish that there are finite-size 
corrections to the contact number and moduli above the 
jamming transition. We also reveal novel finite-size be- 
havior close to the transition that can be scaled to col- 
lapse onto a single curve, firmly establishing the connec- 
tion between jamming and phase transitions. While pre- 
vious work by Olsson and Teitel [8] demonstrated scaling 
collapse in the unjammed regime, our focus is on jammed 
systems above the transition. We find that all finite-size 
effects scale with L d ~ N in d dimensions, where d = 2, 3, 
L is the linear size of the system, and N is the total num- 
ber of particles. Such scaling is expected of a system at 
or above its critical dimension [9] and implies that the 
jamming transition has an upper critical dimension of 2. 
This is consistent with the observation that the power- 
law exponents are the same in 2 and 3 dimensions [4], as 
well as an argument that fluctuations should be unim- 
portant for d > 2 [10]. 

We consider disordered packings of TV frictionless 
spheres at temperature T = and pressure, p, with a 
finite-range, repulsive potential between particles i and 

r- 



energy scale. The effective spring constant between con- 
tacts is fc ff = ^ 9 - ^> [2]. Each packing is relaxed 

to a local energy minimum. We then remove the small 
fraction of "rattler" particles that do not contribute to 
the rigidity of the system [4]. 

Before counting constraints for finite systems, we must 
specify what it means for a system to be jammed. One 
possible definition is that, in the absence of rattlers, 
the only zero-frequency vibrational modes are associated 
with global translation of the system [4] . For N spheres 
in d dimensions, there are dN degrees of freedom and d 
global translations so that dN — d of the degrees of free- 
dom must be constrained. This requires that the number 
of conacts, N coniact > dN — d. Since the contact number 
Z = 2Jv y act | we htain 
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This is the contact number required for the system to 
have no soft modes beyond those corresponding to global 
translations. In the infinite-size limit, this reduces to the 
isostatic condition, Z?° Q — 2d, consistent with previous 
results [3, 4]. However, this definition relics on the choice 
of relevant degrees of freedom. Rattlers, for example, 
have no effect on the elastic properties of a packing but 
contribute d zero modes each if not removed. Similarly, a 
sphere can rotate about its center without any effect on 
the packing. Thus, this definition can break down, as it 
does when generalized to packings of ellipsoids [11-13]. 

A more physical requirement is that system have a pos- 
itive bulk modulus, B. The minimum number of contacts 
needed for a packing of N spheres to have a positive bulk 
modulus, A c ™ n tact , is 
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so that the minimum contact number is: 



only if rij < aij, where is the center-to-center dis- 
tance, an is the sum of their radii, and e = 1 sets the 
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In principle, packings with B > are not forbidden 
from having complicated soft modes. For sphere packings 
(with rattlers removed) we have never observed such non- 
trivial soft modes and therefore assume in the following 
argument that they do not exist for generic packings. In 
this case, at least one additional contact above dN — d is 
required for the system to have a positive bulk modulus. 

The origin of this extra contact can be understood 
by treating the size of the periodic box as a degree of 
freedom. When Z < Z^ Q , there are at most as many 
constraint equations as particle degrees of freedom. If 
there are no nontrivial soft modes, it is possible to sat- 
isfy the constraints nj — o~ij for every contact. Thus, 
by Eq. (1), the total energy and pressure must be zero. 
Since any deformation in the linear regime does not form 
any new contacts, the energy remains constant and the 
bulk modulus, B, must be zero. Therefore at least one 
additional contact is needed for the system to have a pos- 
itive bulk modulus or pressure. This additional contact 
corresponds to the last term in Eq. (3) and leads to Zg 
in Eq. (4). For any positive pressure, the contact number 
should satisfy Z > Zg , and we expect that 
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A third possible definition of jamming is that the 
system have a positive shear modulus, G, as well as 
bulk modulus. Dagois-Bohy, et al. [14] have recently 
shown that packings can be constructed to have pos- 
itive bulk and shear moduli by allowing the shape of 
the box to vary during minimization. In two dimen- 
sions, this introduces 2 extra degrees of freedom for the 
square box to distort to a rhombus or rectangle. In 
d dimensions, there are ^d(d + 1) — 1 degrees of free- 
dom corresponding to the shape of the box. Therefore, 
the extension of our counting argument to such "shear- 
stabilized" packings predicts a minimum contact number 
of Z% G ee 2d - 2d/N + d{d + 1)/N. This exactly agrees 
with the findings of Dagois-Bohy, et al. [14]. 

To test the prediction in Eq. (5) and examine finite- 
size effects, we generated packings of particles with har- 
monic repulsions given by Vfcj) with a = 2 for systems 
ranging from N = 64 to N — 4096. For this potential, 
k c g is independent of (and therefore compression) as 
long as the particles overlap. The relative radii in 2 di- 
mensions were chosen from a flat distribution between 
a ee 1 and 1.4 <r, while in 3 dimensions a bidisperse mix- 
ture of ratio 1:1.4 was used. We fixed the box shape and 
used pressure as the control variable to produce pack- 
ings with a positive bulk modulus. These packings cor- 
respond to what Dagois-Bohy, et al. [14] refer to as the 
"compression-only" ensemble. 

Mechanically stable configurations were generated for 
a range of pressures spanning 7 orders of magnitude. In a 
square (cubic) periodic box, particles were placed at ran- 
dom. The system was then quenched to a local energy 
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FIG. 1. (a) Z — and (b) G as a function of pressure for 
different system sizes in 2 dimensions. For both quantities, 
the power law exponent of 1/2, observed at high pressures, 
agrees with the known scaling for harmonic potentials. At low 
pressures, however, finite-size effects dominate. G is averaged 
over configurations and shear directions. 



minimum (using a combination of linesearch methods, 
Newton's method and the FIRE algorithm [15] to max- 
imize accuracy and efficiency), and the packing fraction 
was adjusted until a target pressure was reached. Sys- 
tems were thrown out if the minimization algorithms did 
not converge. For each dimension, system size, and pres- 
sure, quantities were averaged over 1000 to 5000 pack- 
ings. The shear and bulk moduli were calculated via 
linear response from the dynamical matrix as in [16, 17]. 
In finite systems, there is a well-defined linear regime in 
which the contact network is fixed [18]. By using linear 
response, we ensure that the elastic moduli are calculated 
in this regime. 

Fig. 1(a) and (b) shows both Z — Z^ Q and G versus p 
in 2 dimensions. Similar results are obtained in 3 dimen- 
sions (see Figs. 2 and 3). As expected, Z — Z^ a ~ p 1 / 2 at 
high pressures, consistent with previous studies [3, 4], but 
approaches 2/N at low pressures in accord with Eqs. (4) 
and (5). Thus, one extra contact is needed beyond the 
isostatic value in order for the bulk modulus to be posi- 
tive, as predicted. Moreover, Fig. 2(a) shows that the 
data collapse when (Z — Z^ ) N, related to the total 
number of excess contacts, is plotted versus p x / 2 N . 

It is not obvious from constraint-counting arguments 
alone that at the jamming transition the contact num- 
ber data should obey finite-size scaling: Z — Z(£ = 
N y F(p z N) for some y and z. However, if it does then 



FIG. 2. Collapse of (a) Z — Z^ a and (b) G in 2 and 3 dimensions. In the zero pressure limit, (Z - Z? ) TV -> 2 (dashed 
line), which corresponds to a single contact above isostaticity. (c) Collapse of Z — Zg in 2 and 3 dimensions, (d) Collapse of 
G — Go in 2 and 3 dimensions. The scaling function is qualitatively similar to that of Z — Zg . Inset: the plateau value Go is 
proportional to Symbols and colors are the same as in Fig. 1. 



we can show that y = — 1 and z = 1/2, consistent with 
Fig. 2(a). By counting constraints, we have argued that 
Z - Z? ->■ Z% - Z? = 2 /TV as p -> 0. This is satisfied 
if lim 2; _j.o -^X^) = 2 and y = — 1. In the large TV limit, on 
the other hand, we must recover the asymptotic scaling 
relation Z — Z^ a ~ p 1 ^ 2 , independent of TV. This can 
only be satisfied if F(x) ~ x at large x, and z = 1/2. 
Therefore, the only possible scaling is 

Z-Z£ = ±F(pV 2 N), (6) 

where F(x) ~ 1 for small x and -F(x) ~ x for large a: (see 
Fig. 2(a)). 

The shear modulus, shown in Fig. 1(b), displays G ~ 
p 1 / 2 at high pressures, again consistent with previous 
studies [3, 4]. As the pressure is lowered, however, G 
develops a plateau that is proportional to jj. For a sys- 
tem of TV spheres, one would expect that if one extra 
contact is required to constrain the size of the periodic 
box so that B > 0, additional contacts would be required 
to constrain the shape of the box as well so that G > 0, 
as found by Dagois-Bohy, et al. [14]. However, Fig. 1(b) 
shows that although the shear modulus is not positive 
in all directions for all configurations, the angle- and 
configuration-averaged shear modulus is positive with the 
addition of only one extra contact. To understand this, 
note that the shear modulus measures the response to a 
deformation at constant volume; the size of the periodic 
box is held fixed under shear strain and is therefore no 
longer an independent degree of freedom as it was under 
compression. This allows the extra contact in Eq. (3) to 



do double duty- it can contribute to the stability of the 
system against shear as well as compression. This extra 
contact is the origin of the plateau in G. 

Fig. 2(b) shows that, like (Z ~ Z^ ) TV, GN also shows 
finite-size scaling as a function of p x / 2 N for different 
system sizes and pressures. Note that the slight TV- 
dependence for large pressure in Fig. 1(b) completely 
vanishes when the data is scaled (Fig. 2(b)). This is a 
result of the non-trivial scaling function at intermediate 
pV 2 N. 

The plateaus at low p x / 2 N in the scaling functions for 
(Z - Zg Q ) TV and GTV result from the fact that Z? con- 
tacts per particle are not enough for the system to have 
a positive bulk modulus - one additional contact is re- 
quired. These plateaus can be subtracted off in order to 
study the system-size dependence in greater detail. In 
this case, Fig. 2(c) shows that Z — > Zg at low pressures, 
confirming Eq. (5). Importantly, as we asserted above, 
no properly minimized configurations are observed that 
satisfy both Z < Z^ and 73 > 0. 

Note that (Z - Z%) TV, like (Z - Zg Q ) TV, collapses 
onto a single curve when plotted versus p x l 2 N (Fig. 2(c)). 
There is, however, a crossover to Z — Zg ~ pN for 
p 1 / 2 N < 0(1) in both 2 and 3 dimensions. This scal- 
ing arises because quantities like Z — Z^ Q should only 
be singular at the jamming transition at p = in the 
thermodynamic limit; in finite systems they should be 
analytic around p — 0. Given the existence of scaling 
collapse, which has the form of Eq. (6), the first two 
terms in the Taylor expansion of Z — Z^ Q in powers of p 
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FIG. 3. GN ~ [Z - Zg„) N in both 2 and 3 dimensions. 
Symbols and colors are the same as in Fig. 1. 

must be 

Z-Z&a^+dpN, (7) 

with constants cq and c\. This is precisely what we ob- 
serve in Fig. 2(c), where cq — 2 reflects the extra contact 
at the transition. 

For the same reason, we find the same crossover in 
the shear modulus when we subtract the plateau value, 
G ~ ^. Fig. 2(d) shows that (G - G ) N again col- 
lapses in both 2 and 3 dimensions when plotted against 
p 1 / 2 N, with G - G ~ pN for p x / 2 N < 0(1). The qual- 
itative similarity between {Z — ZjJ) N and (G — Go) N 
underscores the dependence of the shear modulus on the 
contact number. Indeed, we find that to a very good 
approximation, GN is a pure power law in (Z — Z£ ) N 
(Fig. 3) , in accord with recent results of Dagois-Bohy, et 
al. [14]. 

We have also studied the finite-size scaling of the bulk 
modulus, B, which scales as p° for harmonic repulsions. 
Therefore, B approaches a constant value, Bo, as p — » 0. 
As with the coordination number and shear modulus, we 
subtracted off the plateau value to study B — Bo. The 
quantity B—Bq is necessarily quite sensitive to Bo, which 
is large, in contrast to Go, which is of order 1/N. Our 
results are consistent with (B — Bo) N collapsing in both 
2 and 3 dimensions as a function of p x / 2 N with the same 
asymptotic behavior as {Z — Zg) N and (G — Go) N. 

Discussion. We have argued that an appropriate def- 
inition of jamming is that a system can support an ex- 
ternal stress. One could either restrict the definition to 
a compressive stress, requiring B > 0, or to any stress, 
requiring B > and G > 0. If one requires B > 0, 
then sphere packings require one additional contact in 
the entire system, beyond the number calculated for the 
isostatic condition, in order to become jammed. If one re- 
quires both B > and G > 0, then d(d + l)/2 additional 
contacts are required. 



Our results provide a simple interpretation of the re- 
sults of Moukarzel [19], who found that the elastic mod- 
uli vanish in the large N limit for random networks with 
Z = 4 in d = 2. Comparing Z = 4 to Z? Q (Eq. 2), we see 
that Z > Zg so that Z-Z£ Q = A/N. For random spring 
networks, the bulk and shear moduli scale with Z — Z^ D , 
implying that B and G both scale as 1/N for Z = A. 
Thus, our constraint counting arguments imply that B 
and G should vanish as 1/N as N — > oo, consistent with 
Moukarzel's results. 

We find that Z - Zg 0) Z - Z% and G are analytic 
around the jamming transition in finite systems and ex- 
hibit finite-size scaling collapse, a defining characteristic 
of phase transitions. These results cannot be understood 
from constraint counting alone, and provide direct evi- 
dence that the jamming transition is a phase transition. 

The finite-size scaling that we observe depends on the 
total number of particles, TV, rather than on the sys- 
tem length, L, in both 2 and 3 dimensions. For first- 
order transitions, quantities exhibit scaling collapse with 
N ~ L d 1 the number of particles in the system, not with 
L, the linear size of the system [21]. For second-order 
transitions in systems at or above the upper critical di- 
mension, finite-size scaling also leads to collapse with 
N [9, 22]. Earlier observations that the exponents do 
not depend on dimension in d = 2 and 3 [2] and an Imry- 
Ma-type argument of Wyart [20] both suggest that the 
jamming transition has an upper critical dimension of 2. 
Our result that quantities exhibit scaling collapse as a 
function of p 1 / 2 N is therefore consistent with both the 
first- and mean-field second-order character of the jam- 
ming transition. 
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